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Abstract. We consider a model of heat conduction networks consisting of oscillators in 
contact with heat baths at different temperatures. Our aim is to generalize the results 
concerning the existence and uniqueness of the stationnary state already obtained when the 
network is reduced to a chain of particles. Using Lasalle's principle, we establish a condition 
on the disposition of the heat baths among the network that ensures the uniqueness of the 
invariant measure. We will show that this condition is sharp when the oscillators are linear. 
Moreover, when the interaction between the particles is stronger than the pinning, we prove 
that this condition implies the existence of the invariant measure. 

1. Definitions and Results 

1.1. The motivations. We consider an arbitrary graph. At each vertex of this graph, there 
is a particle interacting with the substrate and with its neighbours. Among these particles, 
some are linked to heat baths; an Ornstein-Uhlcnbeck process models this interaction. Given 
this graph, we establish conditions on the disposition of the heat baths that entails existence 
and uniqueness of the invariant measure. 

When the graph is reduced to a chain, each extremal particle is connected to a heat bath. 
This model has been studied in [EPRB99b, EPRB99a, EHOO, RBT02]. The uniqueness of 
the invariant measure is obtained using controllability properties. This property is deeply 
connected to the geometry of the chain: the behaviour of the extremal particle entails the 
behaviour of its neighbour and so on. . . The existence of the invariant measure when the 
interaction is stronger than the pinning has also been obtained. These results were used 
in [Car07] to solve some variations of this model developed in [BOOS] on the one side and [LS04] 
on the other. 

To avoid the particular geometry of the chain, we work with general networks. These 
heat conduction networks have been introduced in [MNV03] and [RB03]. Let us notice that 
an Ornstein-Uhlenbcck process is the sum of a damping term and an excitation term. To 
understand the effect of each of these quantities, we will not suppose that the heat baths have 
non-negative temperatures. A recent work of [BLLO08] uses this kind of result to prove the 
existence of a self-consistent temperature profile. We will see that their results are closely 
related to the geometry of the network they consider. First we introduce the model. Then we 
will state our main results linking existence and uniqueness of the invariant measure to the 
disposition of the heat baths. Intuitively, the existence and the uniqueness of the invariant 
measure is related to the disposition of the damped particles, i.e. the particles interacting with 
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a heat bath. The disposition of the heat baths at non-negative temperatures has an effect on 
the regularity of the invariant measure. 

1.2. The model. Let G ~ {V,^) be a connected graph with vertex set V. Two vertices 
i, j (z V are nearest neighbours if there is an edge between them: i ^ j. Every node i e V 
holds a particle of unit mass described by its position qi and momentum pi. Each particle 
is pinned and interacts with its neighbours. The pinning potential V and the interaction 
potential U are both assumed to be convex polynomial functions. The Hamiltonian of this 
system is the sum of kinetic and potential energies 

H{q,p) = Y.\lp^^+ ^fe) + ^ E ^(^^ - ■ 

Among the particles, some of them are damped; the non-empty subset I? C V is called 
the damped set. Among the damped particles, there is a non-empty subset dV C V called 
the boundary set. Each particle i G dV of the boundary set is connected to a heat bath at 
(non-negative) temperature Tj. The excitation transmitted by the heat bath is modeled by 
a Brownian motion. Thus, the dynamics of the heat conduction network is described by the 
Hamiltonian system of stochastic differential equations, for any z G V, 

, . f dq.^ = Pi dt, 

^' [dp., = -dq^H dt - pa.ev dt + ^/Wll.^av dB,, 

where {-Bi, i G 9V} are independent Brownian motions. 

Remark. Damped particles can be considered as connected to heat baths at possibly null 
temperatures. 

In the following we write N ~ \V\ the number of particles and n ~ 2N. We study the 
diffusion Z = {q,p), solution of the stochastic differential equation (1) via its semigroup (Pt )- 
The adjoint of Pt acts on probability measures: for any measurable set A and probability 
measure /j,, 

P,V(v4)= / Pt(z,A)/i(dz). 

Finally, £ denotes the generator of the diffusion (1). Recall that for any smooth function / in 
the domain Vc of £, 

= {H, /} - J2p.dpJ+ J2 Trdlf, 
iev iedv 

where {•, •} denotes Poisson bracket. The operator C* is the formal adjoint of C. 

We are interested in the existence and uniqueness of a stationary state for heat conduction 
networks, i.e. in probability measures /i such that for any t > 0, Pt*/i = /i, or equivalently, 
= 0. 

1.3. The main results. First of all we will recall and complete some results on the harmonic 
case, when U{q) = V{q) ~ 9^/2. In this setting, the stochastic differential equation (1) 
describing the system is linear. 

In the following we will denote (ei)ig{i_...^„} the canonical base of E". A natural space 
occurring in heat conduction networks is the controllability space (see e.g. [Won79]). Let M 
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Figure 1.1. Example of a heat conduction network : Damped particles are 
striped, boundary particles are drawn in black 

be an ri X n matrix and X C V a subset of particles. We will denote £m,i the smallest vector 
space spanned by (ei+Ar, i e X) stable by M, i.e. 

£m,-i = Span {Af'^Cj+AT, iel,k&n]. 

We say that the graph is asymmetric if dim f a/.p = n, where M will be defined in equa- 
tion (2). Examples suggest that this condition is related to an asymmetric disposition of the 
damped particles in the network. Our first theorem states that the existence and uniqueness 
of the invariant measure is determined by the position of the damped particles. 

Theorem 1. When the potentials are harmonic, there exists a unique invariant measure if 
and only if the network is asymmetric, i.e. dim Sm^v = n. 

When this condition is not satisfied, i.e. dim Sm^v n, there exists a quantity invariant 
with respect to the Hamiltonian flow and the invariant measure is not unique. 

Let us notice that this uniqueness condition was already known (see [EZ04] ) . We give here 
a new proof using completeness and provide an explicit quantity invariant by the Hamiltonian 
flow. 

The shape of the support of the invariant measure is then described by the position of the 
heat baths. 

Theorem 2. Assume there exists a unique invariant measure. The invariant measure has full 
support if and only i/ dim £]\i,dv ~ 

Then we will be interested in the anharmonic case where potentials are only supposed to be 
convex polynomials. We will first state a theorem where all the damped particles are excited, 
i.e. V = dV. Hormander's condition will be defined precisely in Section 3.1. The following 
theorem will be a straightforward consequence of the weak controllability result obtained 
in [Hai05a]. 

Theorem 3. If V ~ dV and Hormander's condition is satisfied, there exists at most one 
invariant measure. 

Finally, in the more general setting where dV C V, wc will provide a condition (see Sec- 
tion 3.2. Condition (5)) on the disposition of the damped particles that entails the uniqueness 
of the invariant measure. We will show that this condition reduces to the asymmetric dispo- 
sition of the heat baths when the potentials are harmonic. This condition will be established 
using a dynamic description of the diffusion through Lasalle's principle. This condition (5) 
states that when the damped particles are fixed, the system can only be in the equilibrium 
position. Recall that H is convex, thus H reaches its minimum at an equilibrium point cq. 
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Theorem 4. When the diffusion is asymptotically strong Feller at cq and the stability condi- 
tion (5) is satisfied, the heat conduction network can have at most one invariant measure. 

Notice that this condition docs not dependent on the temperatures of the heat baths. We 
will discuss the asymptotic strong Feller condition in Appendix A. We will prove that this 
regularity condition is satisfied in the harmonic setting as soon as the graph is asymmetric. 
We could not reach the generality of anharmonic potentials but we expect that this property 
is still true in this setting. 

The importance of the stability condition is strengthened by the following fact about ex- 
istence of the invariant measure. In [RBT02], it is shown that when the pinning potential 
is weaker than the interaction potential and the network is a chain, there exists an invariant 
measure. To generalize the method to general networks, we will also have to suppose that the 
stability condition is satisfied for a limit Hamiltonian (see Section 4, Condition (6)). We will 
call this condition the rigidity condition. Let u (rcsp. v) be the degree of the polynomial U 
(resp. V). 

Theorem 5. If the rigidity condition (6) is satisfied and u> v, i.e. the interaction is stronger 
than the pinning, then there exists a unique invariant measure. 

Section 2 is devoted to the proofs of Theorems 1 and 2. In Section 3, we will briefly prove 
Theorem 3 using results obtained by M. Hairer in [IIai05b]. Then, we will present Lasalle's 
principle and the associated stability condition. Finally, we will present the proof of the 
existence of the invariant measure when interaction is stronger than pinning in Section 4. 



2. The harmonic case 

When the potentials U, V are harmonic, the system (1) of stochastic differential equations 
is linear. We describe the network using the adjacency matrix A = {Sin^j)i_j and the degree 
matrix D ~ diag(^j-^j 1). The Laplace operator A on the graph G is A = Z) — A. We denote 

T = I + A, Ix> = diag(ligx') and Tgy — diag(-y/27^1iggv)- Then, the stochastic differential 
equation (1) can be written: 

(2) dZt= MZtdt + adBt, 

with 





( 






V 


" ) 



2.1. Existence and Uniqueness: First part of Theorem 1. Usually, the existence of the 
invariant measure is obtained via a compactness property. We provide a new proof in the 
harmonic case using the fixed point theorem on complete spaces. We could not generalize this 
argument to a more general setting. 

We say that a matrix M is stable if its eigenvalues have strictly negative real part. We 
begin with two linear algebra lemmas proving that when the matrix a is null, the noise-free 
system is contracting if and only if the asymmetry condition dim £m.v — n is satisfied. 

Lemma 2.1 (Projection). Let q e and dim £m,v = n. If for every fc e I-v^^q = 
then q ^ 0. 
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Proof. Let v G £m,v- There exist Xik such that v ~ J2iev k ^ifcT'^ei. Thus, 

ieV.k 

= 

and q E j,. FinaUy, since jy — {0}, we get q — Q. □ 
Lemma 2.2 (StabiUty). //dim £m.v ~ n then M is stable. 

Proof. We have to prove that under the asymmetry condition, the eigenvalues of M have 
strictly negative real part. 

The eigenvector {q,p) ^ associated to the eigenvalue A satisfies 

p = Xq 
-Tq - It>P = Xp. 

Replacing p by Xq and doing scalar product with q, we get 

X^\\qf + X\\Ij,qf + {Tq,q)^Q. 

Thus, the eigenvalues A are either negative or conjugated (with strictly negative real part). 
However, the real part of A is null if and only if \\I-Dq\\ = 0. Thus /pp = 0, and Iv^q = 0. 
To obtain a contradiction, we multiply the first equation by T to obtain Iv^P = 0. Then, by 
induction on fc, Iv^'^q = Iv^'^P = 0. 

Finally, using projection Lemma 2.1, q = Q and p = 0. Since this is impossible, M is 
stable. □ 

We now use a completeness result to prove the existence statement of Theorem 1. Let 
(X, d) be a Polish space, V the set of probability measures on the Borel cr-field of X. For any 
fi, V G V, we define 

(3) Wi^l,,y)= inf / j d{x,y)g[dx,dv), 

where C(/i, i^) is the set of couplings of [Htv). The functional W defines the Wasserstein 
distance on the set Vi of probability measures such that for one zg S X, W{5zoiP) < +oo. 
Let us recall (see e.g. [Bol07]) that {Vi,W) is a complete space. In the following, wc denote 
d{n,iy) = W{^l,v). 

Proof of the first part of Theorem 1. We first prove a contraction inequality for Dirac mea- 
sures. Indeed, using a trivial coupling, there exists a non-negative constant such that for 
any starting points x and y: 

\\Pt5,-P:5y\\ = M f f\\u-v\\g{du,dv) 



< / / \\u-v\\Pt{x,du)(^Ptiy,dv) 



= nzt~zy\\ 

< ||e^*x-e^"y|| 

< e-"°*|la;-y||. 



since stability Lemma 2.2 ensures that the matrix M is stable. 
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We finally use a classical argument (see [HM06b, Theorem 2.5]) to control the distance 
between two distinct probability measures /i, v 

\\p:^,^p:v\\<e-^''%^i-v\\. 

Thus, there exists a € (0, 1) such that 

\\Pt^l-Ptl,\\<a\\^l-v\\. 

The functional ^ i-^ -f**/^ is contracting on a complete space. Thus, there exists a unique 
invariant measure ^ such that 

and /i is the unique invariant measure of diffusion (1). □ 

Example 1. One can easily check that the network described in figure 2.2 satisfies the asym- 
metry condition but does not satisfy the condition described in [MNV03, Section 2.3, 2.4]. The 
damped particles are drawn in black. 



1 




4 



Figure 2.2. An example of an asymmetric network with 5 particles 
2.2. The non-uniqueness: Second Part of Theorem 1. 

Lemma 2.3. // the graph is not asymmetric, i.e. dim £m^t> n, then there exists a quantity 
K invariant by the Hamiltonian flow and independent of the behaviour of the particles on the 
boundary set dV. 

Proof. We look for a quantity K such that 

K{q,p) = a{z,q)^ + {z,pf, 

where (z, q) ~ X^iev denotes the usual scalar product, a is a real and z a vector. We 
look for a, z such that K is independent of {pi, i € T)) to obtain CtK — and such that 
{H, K} = 0. Since F is symmetric, we easily obtain 

{H,K}^2{z,p){az~Tz,q). 

Since dim (fj^p) — 1; T has an eigenvector z e S^jxi with eigenvalue a e M. Then K{q,p) is 
independent of pi, i € V. □ 

Proof of the second part of Theorem 1. For any smooth function / in , wc decompose the 
formal adjoint of £, 

^ lev i€dV 

-{HJ}+CTf. 
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Let fi be an invariant measure. Using Lemma 2.3, there exists a polynomial function K such 
that for any 7 > 0, fi-y is invariant where 

and Z-y is the normalizing constant. □ 

Example 2. When the network is the diamond described in figure 2.3 (the damped particles 
are drawn in black), we obtain the counter-example given in [MNV03, Section 2.3]. Indeed, an 
eigenvector is (0,1,0,-1) with eigenvalue 3. Thus, the quantity 

K{q,p) = 3 + 

is invariant by the Hamiltonian flow and independent 0/ (qi, ^3,^1,^3). 



1 
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Figure 2.3. An example of a non- asymmetric heat conduction network 

2.3. Support of the invariant measure: Proof of Theorem 2. To complete the study 
of the invariant measure, we propose a characterization of its support. Since the pinning and 
coupling potentials are harmonic, let us recall that the diffusion process starting from 2 G 
can be written at time t > by 

Z^=e^''z+ f e^'^'-'^adBs. 
Jo 

Thus, the covariance Kf of the Gaussian process is Kf = e^^'^*~'*V(T*e*^*'^*~*-' ds, 

where M* denotes the adjoint of M. Thus, {Kt)t satisfies the differential equation 

dtKf = aa* + MKf + KtM* . 

Proof of Theorem 2. Let us suppose that dim Em.ti = n. Then, uniqueness Theorem 1 states 
that the invariant measure is unique. Let us consider the nxn matrix Q solution of Lyapunov 
equation 

(4) MQ + QM* = -a(j\ 

Then, if Kq ~ Q, for any t > 0, Kt ~ Q and Q is the covariance of the invariant measure. 
Then, the matrix Q can be written as 

Q= r e^''aa*e^''*dt. 
Jo 



8 



ALAIN CAMANES 



Moreover, using Lemma 2.3 in [SZ70], 

rank(Q) = dim £M,dv- 

That finishes the proof of Theorem 2. □ 

Let us notice that the uniqueness of this solution using Theorem 2.2 in [SZ70] gives an 
alternative proof of uniqueness Theorem 1. Moreover, in [EZ04], Lyapunov equation (4) is 
used to compute the direction of the heat current in simple networks. 

Remark. If dim £M,dv 7^ then the semigroup is not strong Feller. Indeed, the linearity of 
this equation implies that, for any z G Sj^ gy, G Sj^i gy. Thus, Pt{z^£-^ ^y) = 1^^ ^^{z) is 
not continuous. This remark has to be linked with the full support property of the invariant 
measure in Theorem 3. 

2.4. Example. To end this section, we present a heat conduction network where all the 
particles are damped and the unique invariant measure has a degenerate support. 

In the case described by figure 2.4, the network is made of 6 particles. Tedious but straight- 
forward computations give dim Em.v ~ 12 but dim Smmv = 10. The damped particles are 
represented by striped circles, the boundary set by black ones. 



2 




Figure 2.4. An example of a heat conduction network with a unique invari- 
ant measure with degenerate support 



Remark. In view of [BLLO08] , it could be interesting to wonder if there exists a heat conduction 
network where the solution Q to the Lyapunov equation satisfies Qi+N, i+N = for some i e V. 
In this case, under the invariant measure, particle i would be motionless. We couldn't find 
such a graph. 

3. Uniqueness of the invariant measure 

As we could see in the previous section, when the potentials are harmonic, the linearity al- 
lows to do explicit computations and we can describe precisely the invariant measures. When 
potentials are anharmonic, such a description is no more possible. However, we can generalize 
uniqueness results in the following way. Theorem 3 relies on controllability and regularity 
results established via Hormander's condition. Theorem 4 relies on Lasalle's principle. The 
regularity hypothesis are weaker and the temperature conditions are relaxed. As in the har- 
monic setting, the uniqueness of the invariant measure is related to the damped particles 
only. 

3.1. Hormander's condition and uniqueness. In this subsection, we summarize what can 
be done using control. We use the strategy developed in [HaiOSa] to prove the uniqueness of 
the invariant measure when all the damped particles are linked to a heat bath, i.e. V = dV. 
We recall the following regularity property of semigroups and differential operators. 
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Let £* be the formal adjoint of the generator £. Let r be the cardinahty of the boundary 
set dV, ao E M. and Xq, Xi, . . . , Xr be the vector fields such that 

r 

= ao+Xo + Y.^l 

i=l 

We define inductively the Lie algebra Z{z) at point z E K" as the algebra generated by 
{Xi}i=i^,,,^r , {[^i, r , {[Xi,[Xj,Xi^]]}ijk^Q^^r, ■■• 

We say that Hormander's condition is satisfied if for every z £ R" the lie algebra £(z) has full 
rank, i.e. dim 2{z) = n. 

If Hormander's condition is satisfied then the semigroup and the invariant measures have 
smooth densities with respect to Lebesgue measure (see e.g. [RB06, Corollary 7.2]). 

Remark. When the potentials are harmonic, Hormander's condition is equivalent to an asym- 
metric disposition of the heat baths. More precisely, for any z £ M", dim Em^dv ~ dim £]\j,-d = 
dim £,{z) = n. 

Lemma 3.1. When all the temperatures are equal to 1, the Gibbs measure is a full support 
invariant measure, where 

^iH{dz) = ^e-"'^'''> dz 

and Z is the normalizing constant. 

Remark. To obtain C* hh = in this lemma, it is essential that all the damped particles are 
excited. 

Theorem 3 is then a straightforward consequence of the following theorem. 

Theorem 6 (see [Hai05b], Theorem 4.1). If Hormander's condition is satisfied and the Hamil- 
tonian H has compact level sets then the diffusion can have at most one invariant probability 
measure. 

Remark. It is easy to see from M. Hairer's proof and using Stroock-Varadhan's Support The- 
orem that the invariant measure has full support. 

3.2. Lasalle's principle: Uniqueness of the invariant measure. In this section we relate 
the uniqueness of the invariant measure to the disposition of the damped particles. When the 
potentials are harmonic this condition is optimal. Contrary to Section 3.1 we do not suppose 
that all the damped particles are excited, dV C T). 

This method is based on the contraction properties of the noise-free dynamics used in the 
lectures of J. Mattingly [Mat07] to prove uniqueness of the invariant measure for Stochastic 
Navier-Stokes equations. Our aim is to prove that the contraction point is in the support 
of every invariant ergodic measure. Then, since ergodic measures have disjoint supports, we 
prove that the invariant measure is unique. To control the deterministic diffusion, we introduce 
Lasalle's principle. 

We would like to notice that in this section we just suppose that the semigroup is asymp- 
totically strong Feller at the equilibrium point cq (see Appendix A for more details). 

3.2.1. Lasalle's principle and Stroock-Varadhan's theorem. In this section we recall Lasalle's 
principle. This principle is a generalization of Lyapunov's method. When the derivative of the 
Lyapunov function is not definite negative, this principle says that the solutions of a differential 
system have an attractive point (see [Sas99, p. 198]). In the sequel, we write H{zt) = dtH(zt). 
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Definition 3.1 (Invariant set). A set A C M" is called invariant if all the trajectories starting 
from A stay in A, i.e. for any zq G A, for all t > to, 

e A. 

Theorem 7 (Lasallc's principle, see [Sas99], Proposition 5.22). Suppose there exists a function 
H : K" M_|- of elass satisfying the following conditions: for all a > 0, 

(1) 0(j = {z: H{z) < a} is bounded, 

(2) < 0- 
We denote 

S^^ze ila; H{z) = 0| 
and we consider the biggest invariant subset A of S . Then, for any zq € D,a, 

z? ' A- 

t — >OG 

Remarks. * When S is reduced to an equilibrium point and Qa T this principle implies 

that for any starting point, the trajectory converges to the equilibrium point. 

* As the function H is non-positive on O^, fla is invariant. 

* In the Stochastic Navier-Stokes equations (see [Mat07]), one can compute the speed of 
decrease to the equilibrium point of the deterministic system via a Gronwall inequality. 
In the heat conduction networks we couldn't establish such a result. 

Finally, we recall Stroock-Varadhan's support Theorem. We consider a stochastic system 

dZt^F{Zt) dt + a{Zt)o dBt, 

where F, a are smooth functions and o denotes the Stratonovitch integral. For every zq G 
M", to > 0, let 

Sto.zo = l^^to; ^V' e C", zt = zo + J F{zs) ds + J (t{zs)'iIj{s) ds 

where C~ stands for the set of piecewise continuous functions from [0, oo) to M". 

We recall that the support Supp /i of a measure ^ is the set of points z e M" such that for any 

ball B{zo,e) of radius e > centred in zq, fi{B{zo,e)) > 0. 

Theorem 8 (Support Theorem, see [SV72], Section 5). Using the previous notations, 



SuppPto(zo, •) = Sto,zo- 

3.2.2. Uniqueness of the invariant measure. Recall that the Hamiltonian H has a unique 
minimizer cq. By adding a constant term, we can assume that H{co) ~ 0. Diffusion (1) is said 
to satisfy stability condition if the deterministic system 

{(It = Pt 
p, ^ -dg^H 
Pihev = 

has a unique solution given by zt = cq. 

In our way to prove Theorem 4, we are going to show the following theorem. 

Theorem 9. // the stability condition (5) is satisfied then for every invariant measure /i 

Co € Supp//. 
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Since H is continuous at co, we first notice that for all e > there exists 77 > such that 

K,^ := {z; H{z) < 77} C S(co,e). 

Theorem 9 is proved using the following way: from Lasalle's principle, every solution of the 
deterministic system associated to (5) goes to cq. Hence, it hits K,^. As soon as the dynamical 
system enters K^i, it stays there, since K^j is invariant. If one considers all the solutions 
starting from the ball of radius R, we show that after a finite time T all these solutions are in 
Krj. Finally, using Stroock-Varadhan's Support Theorem 8 we prove that for all z in the ball 
S(0, R), Co is in the support of Pt{z, •)• Now, we give a formal proof of these facts. 




Figure 3.5. Dynamic of the noise- free Hamiltonian system 

For any z G M", let be the hitting time of A',, starting from z, i.e. 

=inf{< > 0; z^ e K,^}. 

We denote 

T = sup T^. 

Since the hitting times of an open set are upper semicontinuous, we get 

T < +00. 

Proposition 3.2. With T and Kjj defined as before, for all z e B{0, R), 

Pt{z,K,,)>Q. 

Proof. Straightforward using the support Theorem 8 and the trivial control ip = Q. □ 
Proof of Theorem 9. First notice that we can use Lasalle's principle since 

Thus, the stability condition (5) involves the convergence of the noise-free solutions to the 
equilibrium point. Let fi be an ergodic invariant measure, e > 0. Since fi is non null, there 
exists a ball Br such that ij.{Br) > 0. Since H is continuous, there exists ry such that Kri C 
;B(co,e). Using the previous Proposition 3.2, there exists T such that for all z G Br, 

Pt{z,K^) >0. 
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Then, for any invariant measure fx, 

M(B(co,e)) > fi{K^) 

^ [ Pt (z, K^) fi{dz) 



> 0. 

Finally, we have shown that cq € Supp /i. □ 

We can now prove uniqueness Theorem 4. 

Proof of Theorem 4- Let us recall (see e.g. [DP06]) that the set of invariant measures is the 
convex hull of the set of invariant ergodic measures. Moreover, if (Pj) is asymptotically strong 
Feller at cq, then cq belongs to at most one invariant ergodic measure [HM06a, Theorem 3.16]). 
But, if there are two invariant ergodic measures, the previous Theorem 9 ensures that cq is in 
both supports. That's impossible! 

Finally, there is a unique ergodic invariant measure and, since the convex hull of a point is 
reduced to that point, diffusion (1) has at most one invariant measure. □ 

3.2.3. The harmonic case. To improve our understanding of the stability condition defined 
below, we go back to the harmonic case where the pinning and the interaction potentials are 
quadratic. 

Theorem 10. When the potentials are quadratic, the stability condition is equivalent to the 
asymmetry condition. 

Proof. In this section, we want to prove that if dim Sm^v = n, zt = is the unique solution 
of the system 

qi = Pi 
Pihev = 0. 

Let us rewrite the equation associated to particle i € T>, (z, e^+jv) = 0. But, if we derive this 
equation with respect to the time parameter, as = on T>, 

qj = constant. 

Thus, we can write this equation (recall that A is symmetric) (z, Aei+jv) = 0. By induction, 
we obtain for any k E N, i £ V, t G R+, 

(2;,A'=e,+Ar) = 0. 

Since dim £m,v = n, 
thus, 

(z, Gi+Af) EE 0, Vi e V. 

Finally we obtain qi ~ constant, Vi G {1, . . . , N}. Since is the unique solution of equation 



Span{A'=^ej+Ar, i€V, keN} =R 



N 



H(z) = 0. we obtain 



0. 

□ 
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3.3. Example. The example presented in figure 2.4 is a heat conduction network that satisfies 
the stabihty condition (5) but not Hormander's condition. We suppose that the interaction 
potential U satisfies U{q) = and that the pinning potential is quadratic, V((i) = (f' . Let 
us notice that all the particles are damped. Thus, stability condition is easily satisfied. We 
will see in Section 4 that this system has an invariant measure. However, we are going to see 
below that Hormander's condition is not satisfied. 
First of all, [Xq, i9pj = ~dq^ + dp^ , thus dq^ G £. 
We compute the Lie bracket between Xq and dq-^ , 

[Xo,dq,] = Y^^lq^Hdj,^ 

iev 

= V"{qi)dp, + U"{qi - q3)dp, + U" {qi - q4)dp, 
■■■-U"{q3-qi)dp,-U"{qi-qi)dp,. 

So, 

U"{qi - q3)dp, + U"{qi - qi)dp, E £{z). 
Computing the Lie bracket between the preceding quantity and , we obtain 

U"'{qi - q3)dp, + U"'{qi - q4)dp, G £(z). 

With a similar computation we get, U^^\qi ~ q?,)dp.^ + U^^^qi — q^jdp^ G S,{z), thus 

5p3 + G ^{^-)- 
We finally distinguish the following cases. 

* If qi ^ 93 and qi ^ 54, Vandermonde's determinant does not vanish, hence 

dim £(z) = 12. 

* 11 qi — q-i and q\ ^ 54 then dp^ G S,{z), and the last Lie bracket gives dp^ G £(2), so 

dim £(z) = 12. 

The same computations could be done if we consider the particles 2, 5 or 6. 

* If (71 = (73 = 94 = 95 = qe = q, we only get dp^ + dp^ G £(z). Then, the Lie bracket 
with Xq gives 

9,3 +dq^ G £(z), 

so 

V"{q)dp,+V"{q)dp,eZ. 

Thus, 

dim £(z) = 10. 

Finally, Hormander's condition is not satisfied on the manifold defined by {z G M"; qi = ■ ■ ■ = 

qe}- 

4. Existence of the invariant measure 

Recall that u (resp. v) denotes the degree of the polynomial U (resp. V). Henceforth we 
assume that u > v, i.e. the interaction is stronger than the pinning. This section is devoted 
to the proof of the existence statement of Theorem 5. The proof of this theorem is based 
on the arguments developed in [RBT02]. We give the outlines of the proof to emphasize the 
importance of the rigidity condition. 

To prove the theorem we use the classical Krylov-Bogoliubov method. 

Lemma 4.1. Let us suppose that there exists a Lyapunov function W such that 
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(1) CW{z) < CW{z) for a positive constant C , 

(2) there exists a time to > and a sequence a„ t +oo such that 

hm sup ° . . = 0. 

"^°°{z;lV(z)>a„} W[Z) 

Then the semigroup (Pt) has an invariant measure. 

Proof. Wc divide the proof in two main steps. 

First, let K be the compact set {W^(z) < a„}, with n large enough. Using the hypothesis, 
there exist a G (0, 1), b such that 

Pt„Wiz)<aWiz) + blKiz). 

Thus, using the semigroup property, for any n G N, 

PntoW{z) < a"W{z) + b + ab+--- + a''-'^b + a'^-^blK{z) 

< a"M^(z) + ' ^ 

:W{z) < CW{z), we 

we get 



1 - a 

Since CW{z) < CW{z), we get PtW{z) < e'^^W{z) and choosing n with nto < t < {n + l)to, 



PtWiz) = PntoPt^ntaWiz) 



< e^(*-"*") ( a^'Wiz) + — ^ 

\ 1 — a 

< e'^*" (w{z) 



1 - a 

Thus, supt PtW{z) < +00. 

Second, let z G K". We consider the family of measures defined for any bounded continuous 
function / by 



Jo 



Using the first step, sup^i'tiW) < supf. PtW{z) < +oo. Thus, (ft) is tight and Krylov- 
Bogoliubov argument gives the existence of the invariant measure. □ 

To prove that the conditions of the previous lemma are satisfied in the heat conduction 
network setting, we use the Lyapunov function W{z) ~ e^^'^^\ 

Lemma 4.2. Let /3 G M such that < /? < max(7i, i G (9V)~^. There exists a constant C > 
and b> 1 such that 

w{z) - ' n J ■ 

We recall briefly the proof of this lemma. 
Proof. Let (3 G (^0,max(T,, i G dVy^^. Using Ito's formula, 

H{Z-^) = H{z)+Y,T^t- f Y,Pfds + Mt. 

But (Mt) is a martingale with quadratic variation [M,M]^ ~ 2 l^^i^Q-^Ti-pl ds. Thus, 
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PtWjz) 

el3H{z) 



^PMt-aeZlM..M]t^a!^lM,Mh+l3j:,eovT,t-l3f*j:,^^pi ds 



The inequality is obtained using Holder inequality and exponential martingales. □ 

To show the semigroup property in Krylov-Bogoliubov Lemma 4.1, we use a scaling argu- 
ment to reduce our problem to a noise- free dynamics. 

Let {En) be an energy sequence growing to -foo. We use the scaled Hamiltonian 

Let us consider the following stochastic differential equation 
,„ . ( dqi ^ dp^Hndt 

The solutions of this system converge to the solutions of the Hamiltonian system described 
by the limit Hamiltonian (see [Car07, Section 5]) 

H^{q,p) = E y + ^n=vh\" + \ E(% - q^T■ 

Thus, we consider the noise-free system, for any i e V, 
j-^ ^ / * : = dp Ma 

be rigid if 
(6) 

has a unique solution given by Zt = cq. 

Lemma 4.3. // the system satisfies the rigidity condition, then for any zq G M" such that 
Hoo{zii) ~ 1; for any solution of {Sao) starting from zq, for any t > 0, 

f Y.pids>0. 

Proof. Let us suppose that there exists t > such that 



Pi = - dq^ H ao - hev Pi - 

To conclude, we need the following rigidity condition. The differential system (1) is said to 




■Jo 



Let us recall that we supposed that H is strictly convex and reaches its minimum at cq G M". 
Since the rigidity condition (6) is satisfied, Zt = cq and we get Hao{zo) = i?oo(co) = 0. We 
have thus obtained a contradiction. □ 

The existence Theorem 5 is a consequence of the following lemma. The proof can be found 
in [Car07, Lemma 5.2]. 
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Lemma 4.4. Let Z„ be a solution of the scaled system starting from Zn with H{zn) —> oo. 
Then, there exists a subsequence (^n^.)fe such that for any C > 0, < > 0, 



lim 

k — 'OO 



= 0. 



Remark. Contrary to the results obtained for a chain, we cannot prove using this method that 
the convergence speed to the invariant measure is exponentiaL Indeed, this rate was obtained 
via compactness properties that we could not generalize. However, the results of [HM07] are 
still valid : when the interaction potential is quadratic and the pinning potential is at least 
of degree 4, is in the essential spectrum of the extension of the generator C to the space 
L^{e-^" dpdq), with /3 < min(l/ri, i€dV). 

Appendix A. The Asymptotic Strong Feller property 

In this section we prove that the Asymptotic Strong Feller property is satisfied, if the 
potentials are harmonic, when the position of the damped particles is asymmetric, i.e. 
dim £m_t> ~ n. When the potentials are convex polynomials, we expect that the diffusion (1) 
is Asymptotically Strong Feller at the equilibrium point cq. 

The Asymptotic Strong Feller property was introduced by M. Hairer and J. Mattingly 
in [HM06a] to study the stochastic Navicr-Stokcs equation. We recall briefly the main defini- 
tions. 

Definition A.l (Totally separating system). An increasing sequence {dn) of pseudo- metrics 
is a totally separating system if for all x 7^ y G M", 

dn[x,y) 1. 

Definition A. 2 (Asymptotic Strong Feller). A semigroup [Pt) is called asymptotically strong 
Feller a,t z G K" if there exists a totally separating system of pseudo- metrics (dn) and a 
sequence t„ > such that 

lim limsup sup \\Pt„{z, •) - Pt^iv^ ■)\\d„ = 0. 

7—0 n^oo yeB{z,f) 

Proposition A.l (Gradient property, see [HM06a], Proposition 3.12). Let z G M". Assume 
there exists a non- decreasing function C : R+ — )■ M and two positive sequences "f 00, (5„ | 
such that for all differentiate function cf) with ||0||oo = sup^ I0(-2^)I o,''^^ |j V(/)(z)||tx) finite, for 
any y in a neighbourhood of z, 

|VPt„(/)(2/)| <C(||2/||)(||0||oo + <5„||V0||oo). 

Then the semigroup is asymptotically strong Feller at z. 

A corresponding totally separating system of pseudo-metrics is then 

d„(x,y) = lA^-^. 

On 

Proposition A. 2. When the potentials are harmonic and the graph is asymmetric, i.e. 
dim Sm.v = n, the diffusion semigroup is asymptotically strong Feller at any point in M". 

Proof. We are going to use the stability of the matrix M to prove that this property is true. 
First notice that for any unitary vector ^ G M", 

D{PtC^){z)C = E[Dq^{Zn(] 

= E[iDcj>)iZn-iDZm 

= : E[{D<j>)iZnPt]. 
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If diffusion (Zt) satisfies the stochastic differential equation dZt = f{Zt)dt + adBt, pt 
D(f>{Z^)^ satisfies the differential equation (see [Bas98, p.30]), 

dpt = {Df) {Zf)pt dt, 

Po = 

Thus, for any function (f> with bounded derivatives, we have the upper bound 

\\^PtH^)\\ < ||V0||ooE[||pJ], 

and it remains to show that E [|jpt||] when t — > oo. 

In the heat conduction network setting, we denote the Hessian of the potential energy 



Hess 



Then, we have 



with M{t, z) 



dtpt = M{t,z)pt, 
Pa = 

/ 

-Hess{Z^) ~Iv 

When the potentials are quadratic, Af(i, z) = M is constant and since matrix M is stable 
as soon as the network is asymmetric, we have pt for any ^ if and only if dim Em.v ~ n. 
Thus, if the network is asymmetric, the diffusion is asymptotically strong Feller. □ 

When the potentials are no more harmonic, we have to study the behaviour of a non- 
autonomous linear differential equation. 
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